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Abstract 

The problem of the thermal and magnetic destruction of the critical state in composite superconduc- 
tors is investigated. The initial distributions of temperature and electromagnetic field are assumed 
to be essentially inhomogeneous. The limit of the thermomagnetic instability in quasi-stationary 
approximation is determined. The obtained integral criterion, unlike the analogous criterion for a 
homogeneous temperature profile, is shown to take into account the influence of any part of the 
superconductor on the threshold for critical-state instability. 
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While dealing with instabilities of the critical state in hard superconductors, the charac- 
ter of the temperature distribution T(x, t) and that of the electromagnetic field t) are 
of substantial practical interest [1]. This derives from the fact that thermal and magnetic 
distractions of the critical state caused by Joule self-heating are defined by the initial tem- 
perature and electromagnetic field distributions. Hence, the form of the temperature profile 
may noticeably influence the criteria of critical-state stability with respect to jumps in the 
magnetic flux in a superconductor. Earlier [2], in dealing with this problem, it was usually 
assumed that the spatial distribution of temperature and field were either homogeneous or 
slightly inhomogeneous. However, in reality, physical parameters of superconductors may 
be inhomogeneous along the sample as well as in its cross-sectional plane. Such inhomo- 
geneities can appear due to different physical reasons. First, the vortex structure pinning 
can be inhomogeneous due to the existence of weak bonds in the superconductor. Second, 
inhomogenety of the properties may be caused by their dependence on the magnetic field H . 
Indeed, the field H influences many physical quantities, such as the critical current density 
jc, the differential conductivity cr^, and the heat conductivity k. 

In the present paper, the temperature distribution in the critical state of composite su- 
perconductors is investigated in the quasi-stationary approximation. It is shown that the 
temperature profile can be essentially inhomogeneous, which affects the conditions of initi- 
ation of magnetic flux jumps. 

The evolution of thermal (T) and electromagnetic [E, H) perturbations in superconductors 
is described by a nonlinear heat conduction equation [3], 



a system of Maxwell equations, 



and a critical state equation 
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Here v — i/(T) is the specific heat, k — k{T) is the thermal conductivity respectively; jc is 

— * 

the critical current density and jV is the active current density. 

We use the Bean-London critical state model to describe the jc{T, H) dependence, accord- 
ing to which jc{T) = jo[l — — Tq)] [4], where the parameter a characterizes thermally 
activated weakening of Abrikosov vortex pinning on crystal lattice defects, Jq is the equilib- 
rium current density, and Tq is the temperature of the superconductor. 

The dependence jr{E) in the region of sufficiently strong electric fields {E > Ef] where Ef 
is the limit of the linear region of the current-voltage characteristic of the sample [2]) can 



be approximated by a piecewise-linear function jr ~ CTfE, where (7/ = ^ ~ '^"'^^ is the 

ttHc 

effective conductivity in the flux flow regime and r] is the viscous coefficient, $o = -z — is the 
magnetic flux quantum, (t„ is the conductivity in the normal state, Hc^ is the upper critical 
magnetic held. In the region of the weak flelds {E < Ef), the function jr{E) is nonlinear. 
This nonlinearity is associated with thermally activated creep of the magnetic flux [5] . 

— * 

Let us consider a superconducting sample placed into an external magnetic fleld H — 

—* 

dH 

(0, 0, He) increasing at a constant rate — — = ij=const. According to the Maxwell equation 

at 

(2), a vortex electric field E = (0, E^, 0) is present. Here Hg is the magnitude of the external 
magnetic field and Ee is the magnitude of the back-ground electric field. In accordance with 
the concept of the critical state, the current density and the electric field must be parallel: 

— * — * 

E II j;. The thermal and electromagnetic boundary conditions for the Eqs. (l)-(4) have the 
form 



dx 



dE 
dx 



+ Wo[T{0)-To] = 0, T{L) = To, 

=0 

(5) 

= 0, E{L)^0, 



x=Q 



For the plane geometry and for the boundary conditions H{0) — He, H{L) — 0, the 

cHe 

magnetic field distribution is H{x) — Hg{L — x), where L = - — - is the depth of magnetic 
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flux penetration into the sample and wq is the coefficient of heat transfer to the cooler at 
the equilibrium temperature Tq. 
The condition of apphcability of Eqs. (l)-(4) to the description of the dynamics of evolution 

of thermomagnctic perturbations arc discussed at length in [1]. 

In the quasi-stationary approximation, terms with time derivatives can be neglected in 
Eqs. (l)-(4). This means that the heat transfer from the sample surface compensates the 
energy dissipation arising in the viscous flow of magnetic flux in the medium with an effective 
conductivity cr/. In this approximation, the solution to Eq. (2) has the form 



E^—{L-x). 
c 



(6) 



Upon substituting this expression into Eq. (1) we get an inhomogeneous equation for the 
temperature distribution T{x,t), 



(7) 



Here we introduced the following dimensionless variables 



/(p) = -[l + ra;p] 
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T-Tq 
To 



L — X 



and the dimensionless parameters uj — , and, r = ( — ■. ^ , where r characterizes 

the spatial scale of the temperature profile inhomogeneity in the sample. Solutions to Eq. 
(7) are Airy functions, which can be expressed through Bessel functions of the order 1/3 [6] 

e(p) = c,p"'Ky, i^-p'i') + cpy^h,, (^p^/^) + eo(p), (8) 

eo(p) = p"'Ky, i^-p'l^) j\l + ru;p,]prh/, (^p?/^) dp,- 



p'^%/s {Ip'^') j\l + -pi]pf i^^pT) dp.. 

where Ci and C2 are integration constants, which are determined by the boundary conditions 
to be 



dp 



Ci = 0, C2 



p= — 
r 
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From the Maxwell equation (2), the temperature inhomogeneity parameter can be expressed 
in the form 
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a 



1/3 



(9) 



It is evident that a ~ 1 near the threshold for a flux jump, when 



aHl_ 

47rz/jo 



1, even under the 



quasi-stationary heating condition — ^ << 1; where — is the characteristic time of 



K 



the heat conduction problem. 

Let us estimate the maximum heating temperature in the isothermal case w = — > 1. 

Ij 

The solution to Eq. (7) can be represented in the form 



Q{x) = - po- 

near the point at which the temperature is a maximum, x — Xm [7] . 



(10) 



With solution (10) being approximated near the point x^ — with the help of the 



thermal boundary conditions, the coefficient po can be easily determined to be 
and the temperature can be written as 



L2 



e{x) = e. 



4 / L 



(11) 



Substituting this solution into Eq.(7), the superconductor maximum heating temperature 
due to magnetic flux jumps can be estimated as 



(TfH 



Jo 



H 

1 



[L-x, 



(12) 



For a typical situation, when — << (L — x^) the estimation for 0^ is 

L'^ CK 



Or 



jo + ^^{L - Xm) 

c 



ckTo 



{L — Xm)- 



(13) 



Here, the parameter 7 ~ 1 (for a parabolic temperature profile 7 ~ 8). It is easy to verify 
that for typical values of jo = 10^ A/ cm^, H = lO'* Gs/sek, and L = 0,01 cm the heating is 
sufficiently low: << 1. In the case of poor sample cooling, w — 1 — lOerg/ {cm^sK), the 



©m is 



©. 



0,5; 



CWqTq 

i.e., the heating temperature can be as high as STm = ToQm ~ 2K. One can see that 
in the case of poor sample cooling, the heating can be rather noticeable and infiucnces the 
conditions of the thermomagnetic instability of the critical state in the superconductor. 
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Let us investigate the stability of the critical state with respect to small thermal 5T 
and electromagnetic 5E fluctuations in the quasi-stationary approximation. We represent 
solutions to Eqs. (l)-(4) in the form 



T{x, t) = T{x) + expl^\sT(^) , 



(14) 



E{x,t) = E{x) + exp <^ — 



At" 



where T{x) and E{x) are solutions to the unperturbed equations obtained in the quasi- 
stationary approximation describing the background distributions of temperature and elec- 
tric field in the sample and A is a parameter to be determined. The instability region 
is determined by the condition that ReX > 0. Prom solution (14), one can see that the 

characteristic time of thermal and electromagnetic perturbations tj is of the order of ^k/A. 

ST 5E A 

we obtam the followmg 



Linearizing Eqs. (1)- (4) for small perturbations y 
equations in the quasi-stationary approximation: 



T{x)' E{x) 



« 1 



(15) 



dx'^ cH 



Eliminating the variable 5T between Eqs. (15), we obtain a fourth- order differential 
equation with variable coefficients for the electromagnetic field 5E: 



dHE 



A(1 + t) + 



E.. 



fSE_ 
~d^ 



+ X[Xt - B{z)\5E = Q. 



(16) 



Here, we introduced the following dimensionless variables: 



X 



z — 



K 



j{z), j{z) = afE{z) - a[T{z) - Tq], E{z) 



HL 



1-z), 
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E. 



One should keep in mind that the variable T{z) and E{z) are given by Eq. (8), in which 

p — —{1 — z). Using the relation between 5T and 5E given by Eqs. (15), we write the 
r 

boundary conditions to Eq. (16) in the form: 



dC'SE 
~d^ 



0, 



5E\ 



2 = 1 



2 = 1 



0, 



d^5E 
~d^ 



= -W 



dSE 
dz 



z=0 



0. 



'fSE_ 
~d^ 



XtSE 



2=0 



(17) 



2=0 
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where W = is the dimensionless thermal impedance. 

K 

The condition for the existence of a nontrivial solution to Eq. (16) subject to boundary 
conditions (17) allows one to determine the boundary of the critical-state thermomagnetic 
instability in a superconducting sample. This problem is complicated, and its analytical solu- 
tion cannot be found in a closed form. We will consider the development of thermomagnetic 
instability in the dynamical approximation, which is valid for composite superconductors 
with high value of CfE. 

The dynamical character of the instability development leads to the predominance of heat 

D K 
diffusion over magnetic flux diffusion in the sample: r = — — >> 1 [1], where Dt — — and 

— are the coefficients of the thermal and magnetic diffusion, respectively. 

In this case, as seen from Eq. (14), the characteristic times tj of temperature and electro- 
magnetic field perturbations have to satisfy the inequalities tj » {X « 1) and tj « t^ 

(At >> 1), where t^ — jz- and tm = are the characteristic times of the thermal and 
magnetic diffusion, respectively. 

As well known that ([3, 4]), the effective magnitude of conductivity cr/ is greater in com- 
posite superconductors than the one in hard superconductors. We can assume that induced 
normal current afE compensates the decreasing of critical current jc{T), caused by increas- 
ing of temperature and obviously prevents magnetic flow penetration into the sample. In this 
case we can neglect the moving of magnetic flow. In other words, in composites the develop- 
ment of thermomagnetic instability is accompanied with "slow" vortices, with characteristic 

time of increasing i,- ~ -(^ >> (A << 1) or >> = = L (Ar >> 1). 

A Ar A 

In the approximation (r >> 1, Ar >> 1, A << 1) Eq. (16) is reduced to a lower order 

differential equation 

d''5E 



+ 



SE^O. (18) 



In the case of r >> 1, the instability threshold depends on the electrodynamic boundary 
conditions at the surface of the sample only shghtly. Therefore, the thermal boundary 
conditions at the boundaries of the current-carrying layer {z — 0, z — 1) can be neglected 
and one can keep only the electrodynamic boundary conditions to Eq. (18), bacause under 
the condition of adiabatic instability development in the composite superconductor, the 
instability threshold depends on the conditions of heat removal on the sample surface only 
weakly [1]. 

Multiplying Eq. (18) by dE and integrating the result with respect to z over the interval 
< 2; < 1, we obtain 
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A = . (19) 



1 5E^dz 



where we use the equahty 

}(f5E^^^ ^^..fdSE\ ff6E\\ }fSE\ 



and the boundary conditions. The right-hand side of Eq. (19) has a minimum at A = Ac: 



1 

^ J B{z)SE^dz 

Ac = . (20) 

J SE'^dz 



Then 



^- 'd6EY , 
— — I dz 

dz 



1 ^ J \ az I 

- J B{z)nldz = '—^ (21) 

° j SE^dz 



where 



2 /■ 2 , 

Up^^— , Ur^dz^l. 



SE'^dz 





Since we do not know the function SE{z), we try, following [8], to obtain an integral es- 
timation of the instability growth increment and the low boundary of its occurrence. The 

HL 

behavior of the integrand in Eq. (21) is basically determined by the factor E = (1 — z), 

which is equal to zero at z = 1 (the other factors change more smoothly). Hence, the 
integrand reaches its maximum at 2; = and the upper estimate for Ac is 

K<^. (22) 

It is evident that Ac >> 1 and AcT <<latr<<l. Numerical evaluation gives Ac ~ 
10-1 ^ 10^3 at r = 10^. 
Equations (21) and (22) enable one to write the instability occurrence criterion in the form 
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This criterion essentially depends on the boundary conditions and the functions j{z), E{z), 
T{z) (see, Ref.[7]). Inequality (23) can be strengthened by means of an evaluation, 



d5E 



dz > 



^ j SE^dz (24) 



dz 

^ ' 
which can be easily verified by expanding the function SE{z) in a Fourier series: 

SE{z) = Amcos — ^— . 

Let us now try to strengthen inequality (24) further. For this purpose, we consider the 
integral 

1 ^1 

/ 9{z){nl - l)dz = J g{z){nl - l)dz + J g{z){nl - l)dz. 

zi 

The last term can be represented in the form 

1 1 
/ 9{z){nl; - '^)dz = {g+ - 9-) j 9{z){n% - l)dz, 



taking intermediate values of the g- in the range z < zi and g > g-'m. the range zi < z <1 
outside the integral. It is evident that [7] 

1 

j g{z){nl-l)dz<Q 



or 

1 1 
j g{z)nldz < J g{z)dz. (25) 



With inequality (25), the instability occurrence criterion can be represented in the form 

1 2 

j B{z)dz>'^T. (26) 



Inequality (26), unlike the analogous criterion for a homogeneous temperature profile, has 
an integral character and takes into account the infiuence of each part of the superconductor 
on the threshold for the superconducting-state instability. If condition (26) is satisfied. 



then small fluctuations of temperature 5T and electric field 5E in the superconductor will 
exponentially increase with time. The most probable result of the development of such an 
instability would be a transition from a critical state to a resistive one. 

In conclusion, we note that the description of the critical state applied in this paper is 
based on the BCS microscopic theory. Then, a continual approximation is used and the 
physical parameters are assumed to vary slowly at distances of the order of the average 
distance d between vortices. The continuity condition is, therefore, L » d. 

One more limitation occurs from possible superconductor overheating above the critical 
temperature T^, where Eqs. (l)-(4) are not valid. This case is realized if the condition 
je(To, E) > a{T, - To) is fulfilled. 
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